Localization in lattice QCD (with emphasis on practical implications) 
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When Anderson localization takes place in a quenched disordered system, a continuous symmetry can be broken 
spontaneously without accompanying Goldstone bosons. Elaborating on this observation we propose a unified, 
microscopic physical picture of the phase diagram of quenched and unquenched QCD with two flavors of Wilson 
fermions. The phase with Goldstone bosons — by definition the Aoki phase — is always identified as the region 
where the mobility edge of the (hermitian) Wilson operator is zero. We then discuss the implications for domain- 
wall and overlap fermions. We conclude that both formulations are valid only well outside the Aoki phase of the 
associated Wilson-operator kernel, because this is where locality and chirality can be both maintained. 



1. Introduction 

Long ago, a conjecture was made for the phase 
diagram of lattice QCD with two flavors of Wil- 
son fermions and a standard plaquette action pQ. 
(For earlier work see ref. This phase dia- 

gram is displayed in Fig. 1. In the Aoki phase (re- 
gion B in the figure), parity and isospin undergo 
spontaneous symmetry breaking (SSB) by a pion 
condensate. 1 According to the Banks-Casher re- 
lation 4 for the case at hand, the pion conden- 
sate is the response to an (infinitesimal) applied 
twisted-mass [115) . 

(tt 3 ) - 2tt P (0) . (1) 

Here p(X) is the spectral density of the hermitian 2 
Wilson operator H(rrio). The twisted-mass term 
is mi i'4>l5 T 3' l P ■ Inside the Aoki phase, the other 
two pions are Goldstone bosons associated with 
SSB of isospin down to a U(l) symmetry. 

*presenter of plenary talk at Lattice 2003, Tsukuba 
1 For go > 0, chiral symmetry is explicit broken by the 
Wilson term, and (iptp') is not an order parameter. In the 
continuum limit chiral symmetry is recovered, and one can 
rotate the pion condensate back to (iptp^, see e.g. ref. HI . 
2 Conventions: H(mo) = 75_D(mo) where the usual Wilson 
operator is D(mo) = D nalve — W — mo. The sum of the 
bare mass mo plus the Wilson term W is a strictly positive 
operator for mo > 0. The super-critical region where zero 
modes may exist is > amo > —8. 



Region C of the phase diagram is interesting 
because, as explained below, this is where one can 
use domain- wall [61 7j and overlap |8l!)j fermions. 
According to the original conjecture Q . all cor- 
relation functions of Wilson fermions are short- 
ranged well inside region C. This was supported 
by quenched numerical results {10. . 3 Today, we 
have convincing evidence that, in region C of 
the quenched theory, there is a non-zero density 
of small-size near-zero modes of H(mo) for any 
go > |12I13| . 4 This raises the following puz- 
zle. Given that p(0) ^ 0, the Banks-Casher rela- 
tion Q implies that the pion condensate is non- 
zero; hence there is SSB of isospin and parity; but 
since isospin is a continuous symmetry, the Gold- 
stone theorem requires the existence of massless 
Goldstone bosons. Thus, the absence of long- 
range correlations appears to contradict the ex- 
istence of a density of near-zero modes in region 
C of the quenched theory! 

Before tackling this puzzle, it is instructive 
to describe the above results using a language 
borrowed from the physics of disordered systems 
(for reviews see ref. J3]). For this purpose we 
may interpret H(mo) as the hamiltonian of a 

3 Exceptional configurations 1111 - defined by the condition 
that H(mo) has a zero mode for some mo - were discarded 
in ref. 1101 . We return to this issue below. 

4 Further evidence comes from the body of numerical work 
using domain- wall and/or overlap fermions. 
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Figure 1. Phase diagram in the bare coupling - 
bare mass plane. The shaded area is the super- 
critical region. 

five-dimensional system, and thus its eigenvalues 
A„ as "energy" eigenvalues. The eigenfunctions 
ty n (x) may be thought of as the wave functions 
of "electrons" that reside on the sites a; of a four- 
dimensional spatial lattice. Disorder is provided 
by the ensemble of gauge-field configurations, in 
terms of which H(mo) is defined. 

Fig. 2 is a cartoon aimed at explaining how the 
small-size, small-A eigenstates of H(mo) emerge 
with increasing disorder. An infinite volume is 
assumed. The first row describes the free op- 
erator Ho(mo), assuming for definitencss that 
amo ~ — 1. The gap is O(l) in lattice units. Next, 
consider one small dislocation (second row), by 
which we mean that the link variables inside a 
small-size hypercube may take any value, but all 
other links are still equal to one. It was shown 
that a bound state - whose eigenvalue is located 
anywhere we like inside the gap - can be produced 
by adjusting the links of the dislocation The 
next step is to consider a dilute gas of small dis- 
locations, choosing the position and shape of the 
individual dislocations at random. A bound state 
produced by a particular dislocation is only neg- 
ligibly affected by all other dislocations. Taken 
together, the (fairly) randomly distributed eigen- 
values of all the bound states fill up the gap. 

The last step is to consider realistic (quenched 
or dynamical) Monte-Carlo configurations. Un- 
like the previous nearly-free cases, now we can- 
not define an "asymptotic value of the potential," 
nor what is the "binding energy." The concept of 
bound states is inadequate and, instead, we have 



exponentially localized states. Likewise, we now 
have extended states instead of scattering states. 
Still, we hypothesize that there are distinct spec- 
tral intervals containing either extended or local- 
ized eigenstates, but not both. 5 

Up to short-distance random fluctuations, the 
mode density of an exponentially localized eigen- 
state behaves like 



l*n(*)l 



exp i 



I 4 



(2) 



On average, the mode density decays exponen- 
tially away from some "center" x^ , where /„ is the 
localization length.® As we increase |A n | starting 
from A„ ~ 0, tunneling becomes easier and, on 
average, l n increases. The localization length di- 
verges when |A n | reaches the mobility edge A c , and 
for | A n | > A c the eigenstates become extended. 7 
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Figure 2. Spectrum of H(mo). On the horizontal 
axis are the eigenvalues X n . The gross features of 
the spectrum are symmetric around A = 0. Part 
of the negative-X spectrum is not shown. 



J We are not aware of any proof of this — widely used — 
assertion. 

6 In four dimensions, typically i S/„(x) will be localized in- 
side a region whose volume is O(Z^). This explains the 
factor of on the right-hand side of eq. Q, for a nor- 
malized eigenstate. 

7 Since J?o( m o) is bounded, we expect another mobility 
edge A c , such that, for A n > A c , the eigenstates are again 
localized (lightly shaded area to the very right in Fig. 2). 
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The value of the mobility edge is a dynamical 
feature. For g <C 1 it will be close to the band 
edge of H (mo), but for larger g it may go down 
and, eventually, even vanish (for artia ~ —1). 

2. Localization, mobility edge and the 
phase diagram 

Let us now reconsider the above puzzle. As- 
sume first that the density of near-zero modes of 
H(mo) arises from extended states only (A c = 0). 
Since these eigenstates spread throughout the en- 
tire lattice, some long-range correlations may be 
expected, and the pion condensate should be ac- 
companied by Goldstone bosons. In contrast, as- 
sume that the condensate arises from exponen- 
tially localized eigenstates only (A c > 0). These 
eigenstates are sensitive only to the link vari- 
ables inside a very small lattice volume (centered 
around and long-range correlations are un- 
likely. We are now closer to resolving the puzzle: 
if p(0) arises from exponentially localized eigen- 
states only, indeed we should not expect any long- 
range correlations. 

What about the Goldstone theorem? The rel- 
evant Ward identity is 

5^(1 + 0(ap)) Pfl f p (p) + 2mif(p) = <7r 3 ) , (3) 

where T(p) and T M (p) are the Fourier transforms 
of (ir + (x) 7r_(y)} and ( J^{x) 7r_ (y)) respectively. 
Jp(x) is a conserved isospin current. (The ± re- 
fer to operators that raise/lower isospin by one 
unit.) As we explain below, when SSB takes 
place in a quenched theory, it is possible to satu- 
rate the Ward identity without Goldstone bosons. 
Instead, the two-point function of the would-be 
Goldstone bosons diverges like 1/mi in the limit 
TOi —* |15I16| . In the unquenched theory, on the 
other hand, the product m\T{p) vanishes in this 
limit for any p ^ 0, and the Goldstone theorem 
applies. 

In the rest of this section we expand on the 
main physical issues. More results, as well 
as missing technical details, may be found in 
ref. PHI, henceforth denoted I. 

After a change of variables ip — 7/175 (see 
also footnote 2), the fermion lagrangian reads 



C = tp (if (mo) — 77771X3)7/;. Let us first consider 
the finite- volume quenched theory, and look for a 
possible 1/mi divergence of the (position space) 
pion two-point function T(x,y). In terms of the 
spectrum of H(mo) we have 



T(a;,y) 



(Afc + 7777i)(A n - imi) 



(4) 



The terms labeled by k (n) refer to the up (down) 
quark. (•} denotes the functional integration over 
the gauge field with the Boltzmann weight of 
the quenched theory. In the denominators, the 
imi terms provide the prescription for integrat- 
ing around the poles at A& = or A„ = for 
mi — > 0. The only terms which may diverge in 
this limit are those where X k = A„, i.e. those 
where k — n (ignoring accidental degeneracies). 
Keeping these terms only we have, for mj —> 0, 

r( X) y) « ±/j2\*n(x)\ 2 \*n(y)\ 2 ^^) 

mi \ ^— ' Ai,+mi/ 



M£|*„(aOI 2 |*„(y)| 2 S(X n ) 



mi 



In a (large but) finite volume, H(mo) has near- 
zero modes practically everywhere inside the 
super-critical region (see I for a more detailed dis- 
cussion). Since the quenched-theory Boltzmann 
weight is strictly positive and independent of mi, 
the pion two-point function indeed has a 1/mi 
divergence! Clearly, (nearly) exceptional config- 
urations 1 1 1 j - the ones responsible for the near- 
zero modes of H(rno) ~ are also directly respon- 
sible for the 1/mi divergence. We conclude that 
the finite-volume, super-critical quenched theory 
is ill-defined for mi = 0, so one must keep mi 7^ 0. 

By summing the previous equation over x and 
then using translation invariance to average over 
y we find, for mi — > 0, 

x * n ' 

which verifies the Ward identity |JS} for p = in 
finite volume. The 1/mi divergence allows for 
p(0) > and, hence, SSB, in the finite-volume 
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quenched theory. 8 

While conceptually important, the previous 
discussion provides little insight into the origin 
of the 1/mi divergence. In fact, we believe that 
in the infinite-volume limit this divergence arises 
only from exponentially localized eigenstates. We 
will illustrate the relevant physics via a heuris- 
tic argument: provided p(0) arises from localized 
states only, there are no Goldstone poles. To this 
end, consider the Fourier transform of the mode 
density H„(p) = a 4 J2 x \^n(x)\ 2 e~ lpx . In view of 
eq. (J2J we assume the ansatz 



n n ( P ) 



1 + pHl 



(6) 



The essential point is that, for p 2 l^ <C 1, one 
should feel only the exponentially decaying enve- 
lope of the mode density (cf. eq. (0) but not the 
short-distance fluctuations. The phase factor in 
the numerator reflects the location of the mode. 
The normalization implies Ti. n (0) = 1. (Positivity 
of the mode density implies \H n (p)\ < 1, so p = 
is the (global) maximum of \Tt n (p)\. This forbids 
a linear term in p in the denominator of eq. ©■) 
For mi — > and p 2 l^ <C 1 we now obtain 



2mir(p) 



V 



27rp(0)(l + o(p 2 I 2 )) , (7) 



where I is the (suitably defined) average localiza- 
tion length of the near-zero modes. Using eq. (0) 

-2 

we conclude that ^^(p) ~ ip^p^) 0(1 ), implying 
that there are no Goldstone poles. 

In conclusion, we arrive at the following phys- 
ical conjecture for toi — > in the super-critical 
region as the simplest one that fits all the exist- 
ing evidence. We define the Aoki phase as the 
region where pseudo-scalar Goldstone bosons ex- 
ist. This region coincides with the part of the 
phase diagram where A c = 0. Inside (and close to 
the boundaries of) the Aoki phase, the chiral la- 
grangian should provide a valid description of the 

8 In the unquenched case one can prove that T(x,y) is 
bounded. The mi — > limit of eq. JS) then yields zero, 
which is the familiar result that there is no SSB in finite 
volume. 



long-range physics |3I17) . Outside the Aoki phase 
A c > 0, and the spectrum of localized eigenstates 
extends down to A = 0. So far, this picture ap- 
plies to both the quenched and the unquenched 
theories. The difference is that, for A c > 0, the 
quenched theory has p(0) > and a divergent 
pion two-point function; the unquenched theory 
has p(A) ~ A 2 for small A due to the fermion 
determinant, and a condensate can only be built 
from extended eigenstates (i.e. when A c = 0). 

3. Domain- wall and overlap fermions 

Domain-wall fermions (DWF) and overlap 
fermions are sophisticated descendants of Wil- 
son fermions, that allow separation of the chi- 
ral and the continuum limits. They are rec- 
ognized today as closely-related realizations of 
the Ginsparg- Wilson (GW) relation JHj- They 
are both built around a Wilson-operator "kernel" 
H(-M), where < a M < 2 is the domain- 
wall height. 9 DWF are five-dimensional Wilson 
fermions, in which only hopping terms in the four 
physical directions couple to four-dimensional 
link variables, which themselves arc independent 
of the fifth coordinate PLTj . The latter takes val- 
ues s = 1,2, . . . ,N S . The left- and right-handed 
components of the quark field are localized on op- 
posite boundaries of the five-dimensional space, 
and chiral symmetry becomes exact in the limit 
N a — > oo [7J. 

Let T(M 1 a$) be the transfer matrix that 
hops the fermions in the fifth dimension. The 
"hamiltonian" H(M,a 5 ) = -\og(T(M,a 5 ))/a 5 is 
closely related to the Wilson operator. One has 
H(M,0) = H(-M). Also, the spectrum of ex- 
act zero modes is independent of a§. We may 
thus study the approach to the chiral limit using 
the physical concepts of the previous section. We 
begin by considering the DWF's PCAC relation 
(the superscript a is an SU(N)-flavor index) [7J 



2m q Jg(x)+2J«(x) 



(8) 



Here A^x) is the partially-conserved DWF axial 
current, d* is the backward derivative, and m q is 
the bare quark mass. The pseudo-scalar density 

9 Here we assume a fifth-dimension lattice spacing as < 1. 
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Jg (x) is localized on the boundaries of the fifth 
dimension and serves as an interpolating field for 
pions. J$ q {x) is another pseudo-scalar density lo- 
cated in the middle of the fifth dimension, which 
gives rise to finite- N s chiral symmetry violations. 
A measure of these violations is the residual mass 
(r, N a ), which may be defined as 

Here we singled out one lattice direction as eu- 
clidean time r. For large r > and large N s , we 
expect, ignoring power corrections, 10 

m res ~ ci exp(-r(f/7 - m n )) + c 2 q Ns . (10) 

This should be true provided A^, > 0, where A'. = 
— log(g)/a 5 is the mobility edge of H(M,a 5 ); 
equivalently, < q < 1. Let us briefly ex- 
plain this result (for the limitations of our anal- 
ysis, see J). The denominator in eq. © is gov- 
erned by a zero-momentum pion and decays like 
exp(— m,T). In the numerator, one should distin- 
guish between the contributions of extended and 
of localized modes of H(M, 05). For large N s , 
modes with eigenvalue near X' c dominate the ex- 
tended modes' contribution. This results in uni- 
versal fifth-dimension wave functions for the left- 
and right-handed quarks, which decay exponen- 
tially with s |19llfi| . The extended modes' contri- 
bution thus factorizes like q Ns exp(-m T r). The 
contribution of localized modes to the numerator 
is not suppressed exponentially with N s , because 
the localized spectrum goes down to zero. This 
contribution can be estimated using the asymp- 
totic form J5J. It will be dominated by localized 
near-zero modes centered around the straight line 
connecting the two space-time points (see / for 
details) and is estimated to behave like exp(— r/I) 
where now I is, in effect, the maximal localization 
length of the near-zero modes. 

According to eq. (|1D|> . m res does not vanish 
exponentially with N s , for any fixed value of 
t. However, let us imagine taking the infinite- 
volume limit while keeping the lattice spacing 

10 In a numerical simulation, the signal of localized modes 
becomes more complicated when full translation invari- 
ance is not enforced. 



fixed. If we extract m res using increasingly large 
r, the first term on the right-hand side of eq. i|lU[l 
will ultimately vanish even for fixed N s , provided 
1/7 > m % . Therefore, the localized modes' con- 
tribution to m res does not necessarily constitute 
a violation of chiral symmetry. 

The overlap operator |8I9| corresponds to the 
double limit N s -> 00 and a 5 -> of DWF. Ex- 
plicitly, 

H(mo) 

aD ov = 1 - 7575 , 75 = — . 11 

\H(m )\ 

This operator satisfies the GW relation |18j . i.e. 
it possesses a modified chiral symmetry (with the 
same algebraic properties as ordinary chiral sym- 
metry). The overlap operator cannot have a fi- 
nite range [20] > and the relevant question is what 
are its localization properties. Exponential local- 
ity of the overlap operator was proved in ref. |21) . 
provided all the plaquette variables are uniformly 
bounded close to one ( "admissibility condition" ) . 
In this case the spectrum of H(mo) has a gap, 
as on the first row of Fig. 2. 11 In realistic MC 
simulations, however, one cannot impose such a 
constraint on the plaquettes. Assuming the mo- 
bility edge satisfies A c > 0, the spectrum looks 
like the last row of Fig. 2. The effect of the lo- 
calized spectrum can be analyzed along similar 
lines. Considering the restriction of % to (say) 
localized modes with |A| < A c /2, denoted 7^, we 
estimate 

(l75 < (^,y)l)<exp(-|x-y|/(20), (12) 

where now I is determined by all eigenmodes with 
|A| < Ac/2. Once the near-zero modes do not 
hamper locality, exponential locality of D ov can 
be established as in ref. |21|. Finally, a similar 
analysis may be applied to the N s — > 00 limit 
of DWF with 05 > 0, which is also described 
by eq. ltTT|l . except that H(m ) is replaced by 
H(M, a 5 ) E2I (see also ref. E3)- 

The main conclusion of the previous discussion 
is the following: DWF and overlap fermions can 
be used only well outside the Aoki phase of their 
Wilson kernel. This situation is quenched- like: 

11 Ref. [211 generalized the proof to the case depicted on 
the second row of Fig. 2. 
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(even) for (dynamical) DWF or overlap simula- 
tions, the Boltzmann weight does not contain the 
determinant of the Wilson operator itself. Nev- 
ertheless, the spectral properties of the Wilson- 
operator kernel are crucial, and, for any ensemble 
of gauge-field configurations, we may (formally) 
introduce two quenched Wilson flavors and look 
for their Aoki phase. 

When the mobility edge of the Wilson kernel is 
0(1) in lattice units, chiral symmetry of DWF 
will be recovered exponentially with increasing 
N s , for large t. 12 Well outside the Aoki phase we 
also expect the average localization length of the 
near-zero modes to be O(l) in lattice units, and 
this guarantees the locality of the overlap opera- 
tor, as well as of the generalized overlap operators 
obtained from DWF with 05 > 0. In contrast, be- 
ing inside the Aoki phase means that the mobility 
edges of both the Wilson kernel and H(M, as) are 
zero. This corresponds to taking the limits q — > 1 
and I — > 00, where eqs. (|10() and (|12|l cease to 
hold. Thus, DWF will develop long-range corre- 
lations in all five dimensions, and the (general- 
ized) overlap operator defined by the N s — > 00 
limit will become non-local, including for 05 — > 0. 

Near the continuum limit, DWF and overlap 
fermions are local. The mobility edge and l/l 
both scale like the lattice cutoff, and p(0) tends to 
zero rapidly for go — > ^2- But, for present-day 
simulations, it is unclear if 1/7 is large compared 
to e.g. the rho-meson mass. How, then, can one 
tell if a given set of simulation parameters lies 
safely outside the Aoki phase? At the very least, 
one should require that l/l will be bigger than 
the mass of all hadrons of interest. The issue 
definitely deserves further study. 

For DWF, eq. I|10|l shows that, by monitor- 
ing the dependence of m res on both N s and r, 
one can extract the crucial spectral properties of 
H{M : a$). The mobility edge is determined in 

12 The situation is more complicated when considering 
the lattice renormalization of the effective Electro- Weak 
hamiltonian. Especially when power-divergent subtrac- 
tions are involved, localized modes may still give rise 
to chiral symmetry violations which are not suppressed 
by any space-time separation. In this case, employing 
smeared links 1241 and/or the "projection method" of 
ref. 1251 may provide better control over the subtractions 
needed to recover chiral symmetry. 



terms of q, and may be extracted from the -/V s de- 
pendence. The (dominant) localization length of 
the near-zero modes may be extracted from the 
t dependence, and the constant c\ provides in- 
formation on the density of the near-zero modes. 
A well-behaved m res (cf. eq. HlOj l') also ensures 
the locality of the effective four-dimensional op- 
erator defined by the N s — > 00 limit. m res is 
routinely calculated in any new DWF simulation. 
For a lattice cutoff a -1 ~ 2 GeV, well-behaved 
m res have been obtained in quenched simulations 
with Iwasaki [251 an( i DBW2 [27] gauge actions. 
Likewise, when overlap fermions are employed, it 
should become a routine practice to determine the 
localization properties of the overlap operator! In 
particular, the relation between the localization 
lengths of the overlap and of the near-zero modes 
should be studied in more detail. 

Dynamical DWF [2B] and overlap simulations 
are, and will be, very expensive. In a dynam- 
ical simulation the danger is that one will end 
up too close to, or even inside, the Aoki phase. 
In this case, one cannot maintain both chirality 
and locality. 13 Having to give up on something, 
we believe that one should insist on locality, at 
the price of doing worse on chirality. The reason 
is that, unlike with approximate chiral symme- 
try (see below), very little is known about how to 
monitor for the physical consequences of deteri- 
orating locality. Any non-locality of an operator 
that satisfies the GW relation should be regarded 
as the source of an unknown systematic error. 

In a situation where the overlap operator (or 
generalizations to 05 > 0) is non-local, locality 
could be maintained, for example, by using DWF 
with modest iV,. 14 For other approximate solu- 
tions of the GW relation see ref. |21j. 15 While vio- 

13 We believe that this applies to present-day thermo- 
dynamical simulations carried out with a lattice cutoff 
a- 1 ~ 1 GeV. 

14 If no pseudo-fermion fields are included, this is strictly 
local. 

15 Local approximations of the overlap exist, too. For ex- 
ample, a rational polynomial approximation may be rep- 
resented as a continued fraction which, in turn, may be 
cast in the form of a five-dimensional action with nearest- 
neighbor coupling 30 . This action is not yet local because 
it contains parameters which are functions of the minimal 
eigenvalue of H(mo), which in turn is a global property of 
the gauge-field configuration. If one fixes these parame- 
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lations of chiral symmetry in Ward identities will 
be non-negligible in this case, the well-established 
renormalization program allows us to subtract 
them. In principle, this is the same situation as 
with ordinary Wilson fermions. In practice, how- 
ever, unlike Wilson fermions, as long 
a few MeVs (or less), it may be possible to carry 
out the subtractions successfully for a wide range 
of weak matrix elements. 
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